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The Coulomb interaction between the two protons is included in the calculation of proton-deuteron 
elastic scattering, radiative proton-deuteron capture and two-body electromagnetic disintegration of 
®He. The hadron dynamics is based on the purely nucleonic charge-dependent (CD) Bonn potential 
and its realistic extension CD Bonn -|- A to a coupled-channel two-baryon potential, allowing for 
single virtual A-isobar excitation. Calculations are done using integral equations in momentum 
space. The screening and renormalization approach is employed for including the Coulomb inter¬ 
action. Convergence of the procedure is found already at moderate screening radii. The reliability 
of the method is demonstrated. The Coulomb effect on observables is seen at low energies for the 
whole kinematic regime. In proton-deuteron elastic scattering at higher energies the Coulomb effect 
is confined to forward scattering angles; the A-isobar effect found previously remains unchanged by 
Coulomb. In electromagnetic reactions Coulomb competes with other effects in a complicated way. 

PACS numbers: 21.30.-x, 21.45.+v, 24.70.-l-s, 25.10.-|-s 


I. INTRODUCTION 

Experimentally, hadronic three-nucleon scattering is 
predominantly studied in proton-deuteron (pd) reactions, 
i.e., in pd elastic scattering and breakup: Proton and 
deuteron beams and targets are available, with and with¬ 
out polarization. The detection of charged particles 
yields complete experiments. In contrast, the charge- 
symmetric neutron-deuteron (nd) reactions are much 
more difficult to perform, since neutron beams are scarce, 
neutron targets non existing, and the detection of two 
neutrons is a complicated experimental endeavor. In 
electromagnetic (e.m.) reactions, proton-deuteron radia¬ 
tive capture has a corresponding advantage over neutron- 
deuteron capture and, furthermore, ^He is a safer target 
with easier detectable breakup products compared with 

3H. 

In contrast, the Coulomb interaction between the two 
protons is a nightmare for the theoretical description of 
three-nucleon reactions. The Coulomb interaction is well 
known, in contrast to the strong two-nucleon and three- 
nucleon potentials mainly studied in three-nucleon scat¬ 
tering. However, due to its I/r behavior, the Coulomb 
interaction does not satisfy the mathematical properties 
required for the formulation of standard scattering the¬ 
ory. When the theoretical description of three-particle 
scattering is attempted in integral form, the Coulomb in¬ 
teraction renders the standard equations ill-defined; the 
kernel of the equations is noncompact. When the theo¬ 
retical description is based on differential equations, the 
asymptotic boundary conditions for the wave function 
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have to be numerically imposed on the trial solutions 
and, in the presence of the Coulomb interaction, those 
boundary conditions are nonstandard. 

There is a long history of theoretical prescriptions for 
the solution of the Coulomb problem in three-particle 
scattering, where different procedures are followed by the 
groups involved. A modified momentum-space integral 
equation approach is used in Refs. nil , whereas the 
configuration-space differential equation approach is used 
in Ref. 0 in a variational framework and in Refs. ii 
in the framework of the Faddeev equations. There are 
more recent formulations lEQ of exact scattering equa¬ 
tions with Coulomb which, however, have not matured 
yet into practical applications. In addition there exist 
approximate schemes: The most brutal one is the de¬ 
scription without Coulomb for the three-nucleon system 
with two protons at those energies and in those kinemat- 
ical regimes in which the Coulomb interaction is believed 
to be irrelevant for observables; such an approximation 
has become standard in recent years Q, and, to our own 
guilt, we admit having used it 0. Reference ^3 extends 
the assumed applicability of that approximation scheme 
by the addition of external Coulomb correction terms to 
those non-Coulomb results. 

In this paper our treatment of the Coulomb interaction 
is based on the ideas proposed in Ref. m for two charged 
particle scattering and extended in Ref. for three- 
particle scattering. The Coulomb potential is screened, 
standard scattering theory for short-range potentials is 
used, and the obtained results are corrected for the un¬ 
screened limit. We rely on Refs. [Hill with respect to 
the mathematical rigor of that procedure. We constrain 
this paper to the description of reactions involving the 
pd system. Thus, we leave out breakup in pd scattering 
and three-body breakup in e.m. reactions with ^He. We 
explain the features of our procedure in order to ease the 
understanding for the uninitiated reader and to point out 
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differences of our treatment relative to Refs. Hi , which 
also are based on Refs. [Hill: 

(1) The calculations of Refs. H Q need improve¬ 
ment with respect to the hadronic interaction. Whereas 
Refs. H 0 limited themselves to the use of low-rank 
separable potentials, we use modern two-nucleon poten¬ 
tials and three-nucleon forces in full without separable 
expansion. In particular, the results of this paper are 
based on the purely nucleonic charge-dependent (CD) 
Bonn potential and on its coupled-channel exten¬ 
sion CD Bonn -|- A 0, allowing for a single virtual A- 
isobar excitation and fitted to the experimental data with 
the same degree of accuracy as CD Bonn itself. In the 
three-nucleon system the A isobar mediates an effective 
three-nucleon force and effective two- and three-nucleon 
currents, both consistent with the underRing effective 
two-nucleon force. A reliable technique Q for solving 
the three-particle Alt-Grassberger-Sandhas (AGS) equa¬ 
tion ^ without Coulomb is at our disposal. We extend 
that technique to include the screened Coulomb potential 
between the protons. Thus, the form of our three-particle 
equations including the screened Coulomb potential is 
completely different from the quasiparticle equations of 
two-body type solved in Refs. 0l3l- 

(2) We work with a Coulomb potential wr, screened 
around the separation r = R between two charged 
baryons. We choose wr in configuration space as 

WR{r) = w{r) (1) 

with the true Coulomb potential w(r) = —, a being the 
fine structure constant and n controlling the smoothness 
of the screening. We prefer to work with a sharper screen¬ 
ing than the Yukawa screening (n = 1) of Refs. 0,13 
We want to ensure that the screened Coulomb poten¬ 
tial Wr approximates well the true Coulomb one w for 
distances r < R and simultaneously vanishes rapidly for 
r > R, providing a comparatively fast convergence of 
the partial-wave expansion. In contrast, the sharp cut¬ 
off (n —> oo) yields an unpleasant oscillatory behavior in 
the momentum-space representation, leading to conver¬ 
gence problems. We find values 3 < n < 6 to provide a 
sufficiently smooth, but at the same time a sufficiently 
rapid screening around r = R; n = 4 is our choice for 
the results of this paper, unless indicated otherwise. The 
screening functions for different n values are compared in 

Fig.n 

(3) The screening radius R is chosen much larger than 
the range of the strong interaction which is of the order 
of the pion wavelength h/rriTrC 1.4 fm. Nevertheless, 
the screened Coulomb potential wr is of short range in 
the sense of scattering theory. Standard scattering the¬ 
ory is therefore applicable. However, the partial-wave 
expansion of the pair interaction requires much higher 
angular momenta than the one of the strong two-nucleon 
potential alone. 

(4) The screening radius R will always remain very 
small compared with the nuclear screening distances 
which are of atomic scale, i.e., 10® fm. Thus, the em- 



FIG. 1: Screening function WR{r)/w{r) as function of the 
proton-proton distance r for characteristic values of the pa¬ 
rameter n in Eq. (Q: n — 1 (dashed-dotted curve) corre¬ 
sponds to Yukawa screening, n = 4 (solid curve) is the choice 
of this paper, and n —> oo (dotted curve) corresponds to a 
sharp cutoff. 

ployed screened Coulomb potential wr is unable to sim¬ 
ulate the physics of nuclear screening properly and even 
more all features of the true Coulomb potential. There¬ 
fore WR is unable to yield the Coulomb scattering ampli¬ 
tude as well as the logarithmic distortion of the Coulomb 
wave function and, consequently, the true Coulomb phase 
shifts. However, since the Coulomb scattering amplitude 
and the Coulomb phase shifts are known and since their 
occurrence in the three-particle scattering amplitudes 
can be spotted, approximate calculations with screened 
Coulomb Wr can be corrected for their shortcomings in 
a controlled way. References 0113 give the prescrip¬ 
tion for the correction procedure which we follow here, 
and that involves the renormalization of the on-shell am¬ 
plitudes in order to get the proper unscreened Coulomb 
limit. 

(5) After the indicated corrections (4), the predictions 
for observables of three-nucleon reactions have to show 
independence from the choice of the screening radius R, 
provided it is chosen sufficiently large. That convergence 
will be our internal criterion for the reliability of our 
Coulomb treatment. 

Section ni describes the practical working of the above 
approximation program in detail. Section EH presents 
some characteristic effects of Coulomb in three-nucleon 
reactions. Section Cvl gives our conclusions. 


II. TREATMENT OF COULOMB 
INTERACTION BETWEEN PROTONS 

Section |H recalled the general idea for including the 
Coulomb interaction in pd scattering and in related e.m. 
reactions by screening and renormalization. This sec¬ 
tion provides the theoretical framework on which we base 
our practical procedure. We are aware that the equa¬ 
tions given here have been develo^d in Ref. 0, but 
their practical realization in Refs. 0,0 differs substan- 
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tially from the work presented here. For completeness 
we rederive some of the equations and explain how we 
solve them; our description aims at elastic pd scattering 
and e.m. reactions involving the pd system. However, 
the essence of the procedure can already be well seen in 
proton-proton (pp) scattering which we therefore use as 
an illustrative example. The example is also an espe¬ 
cially useful test, since exact results for the inclusion of 
the Coulomb interaction are readily available; we recover 
the exact results, easily obtainable in calculations with a 
sharp cutoff Coulomb [l^. 

The numerical results presented in this section refer 
to the coupled-channel potential CD Bonn -|- A which 
allows for single A-isobar excitation. Coulomb acts be¬ 
tween the two protons and between the proton and the 
A"*" in the coupled channel with an isobar. 

A. Proton-proton scattering 

The two protons interact through the strong poten¬ 
tial V and the Coulomb potential w. We introduce the 
full resolvent for the auxiliary situation in which the 
Coulomb potential w is replaced by the screened poten¬ 
tial Wh 

gW(^z) = {z - ho - V - wr)~^, (2a) 

where ho is the kinetic energy operator. The full resol¬ 
vent g^^\z) yields the full scattering state when acting 
on a plane-wave state |pj^) of relative momentum p, en¬ 
ergy e(p), and discrete two-particle quantum numbers v 
and taking the appropriate limit z = e(p) -I- iO. The full 
resolvent therefore also yields the desired S matrix. The 
full resolvent g^^'> (z) depends on the screening radius R 
for Coulomb and that dependence is notationally indi¬ 
cated. Next, we discuss formal manipulations of the full 
resolvent. It can be decomposed according to 

= go{z) + go{z)t^^\z)go{z) (2b) 

with the free resolvent 

go{z) = {z - ho)~'^ (3) 

and the full transition matrix 

t(«)(z) = (u -k wr) -k (u -k WR)go{z)t^^\z). ( 4 ) 

Of course, t^^\z) must contain the pure Coulomb tran¬ 
sition matrix tR^z) derived from the screened Coulomb 
potential alone 

tR{z) =WR + WRgo{z)tR{z), ( 5 ) 

which may be isolated in the full resolvent 

g^^\z) = go{z) + go[z)tR{z)go{z) 

+ - tR{z)\go{z). 


However, an alternative decomposition of the full resol¬ 
vent appears conceptually neater. Instead of correlating 
the plane-wave state |p^) in a single step to the full scat¬ 
tering state by g^^\z), it may be correlated first to a 
screened Coulomb state by the screened Coulomb poten¬ 


tial Wr through 

9 r{z) = {z - ho - wr)~^, (7a) 

9r{z) = 9o{z) + go{z)tR{z)go{z). (7b) 

Thus, the full resolvent can alternatively be decomposed 
into 

9 ^^\z) = gR{z) -k gR{z)f^\z)gR{z), (8a) 

g^^\z) = go{z) + go{z)tR{z)go{z) 

+ 5o(2:){[l + tR{.z)go{,z)\i^^\z) 

x[l +go{z)tR{z)]}go{z) (8b) 

with the short-range operator 

i^^\z)=v + vgR{z)i^^\z). (9) 


Equation Il8bll gives an alternative form for the difference 
of transition matrices [t^^\z) — tfl(z)] in Eq. ®, i.e., 

t^^\z) - tR{z) = [1 -k tR{z)go{z)]i^^\z)[l -k go{z)tR{z)]. 

( 10 ) 

The above equation is the well-known two-potential for¬ 
mula that achieves a clean separation of the full transition 
matrix t^^\z) into a long-range part tfl(z) and a short- 
range part [t^^^{z) — tR{z)]. In this paper the left-hand 
side of Eq. ca is calculated directly from the poten¬ 
tials V and Wr according to Eqs. 0 and ©■ Equation 
m is only introduced by us in order to demonstrate 
that [t^^\z) — tR{z)], even in the infinite R limit, is a 
short-range operator due to the short-range nature of 
V and i^^\z). However, on-shell, it is externally dis¬ 
torted due to the screened Coulomb wave generated by 
[1 + go{z)tR{z)] which together with the long-range part 
tR^z) does not have a proper limit as R ^ oo. This 
difficulty brings about the concept of renormalization of 
on-shell matrix elements of the operators as proposed in 
Refs. PII3 in order to recover the proper results in the 
unscreened Coulomb limit. 

According to Refs. [11113. the pp transition ampli¬ 
tude (pfiyfltlpiRi), referring to the strong potential v and 
the unscreened Coulomb potential w, is obtained via the 
renormalization of the on-shell t^^\z) with z = e{pi) + i0 
in the infinite R limit 

{pfVf\t\p,v,) = lim {Z^^pf){pfiyf\ 

R^oo 

X t^^\e{pi) + iQ)\piVi)Zj^'^{pi)}. 

(11a) 

The transition amplitude {pfvf\t\Pii'i) connects the ini¬ 
tial and final states \piVi) and \pfVf), pf = pi, of the 


4 


considered reaction. However, Eq. (Illall as it stands is 
not suitable for the numerical calculation of the full tran¬ 
sition amplitude; instead, the split of the full transition 
matrix (z) into long- and short-range parts is most 
convenient. For the on-shell screened Coulomb transition 
matrix contained in the limit in Eq. nia 

can be carried out analytically, yielding the true Coulomb 
transition amplitude he., 

{pfiyf\t\p,iyi) = {pflyf\tc\p^lyi) 

+ {Z~^{pf){pfi^f\[t^^\e{pi)+iO) 

R—^oo 

- tR{e{Pi) + iQ)\\piVi)Z^'^{pi)}, 

(11b) 

whereas the limit for the remaining short-range part 
- tR{z)] of the transition matrix t^^\z) has to 
be performed numerically, but it is reached with suffi¬ 
cient accuracy at finite screening radii R. In contrast to 
{pfiyf\tc\PiRi), the short-range part [t^^^z) —can 
be calculated using a partial-wave expansion of Eqs. 0 
and Q. 

The renormalization factor for i? —> oo is a diverging 
phase factor 

Znip) = (12a) 

where (jiRijp)^ though independent of the pp relative or¬ 
bital angular momentum L in the infinite R limit, is given 

by 0 

(I>r(p) = - Vlr{p), (12b) 

with the diverging screened Coulomb phase shift rjLRip) 
corresponding to standard boundary conditions and the 
proper Coulomb one (Jl{p) referring to the logarithmi¬ 
cally distorted proper Coulomb boundary conditions. 
The form of the renormalization phase is readily 

understood by looking back to Eq. m and realizing 
that the external distortion generated by the screened 
Coulomb wave function [l+go{e{p)+iO)tR{e{p)+iO)]\pR) 
carries, in each partial wave, the overall phase fac¬ 
tor ^ 5 . Except for this overall phase factor, 

the screened Coulomb wave approximates well the un¬ 
screened one in the range required by the operator (z) 
in Eq. 113 for distances r < R. Therefore, through the 
renormalization, that unwanted phase factor is changed 
to the appropriate phase factor for the unscreened 

Coulomb wave. 

For the screened Coulomb potential of Eq.JD the in¬ 
finite R limit of 4 > r { p ) is known analytically |lll 

</>ii(p) = «:(p)[ln(2pi?) - C/n], (12c) 

k { p ) = ap/p being the Coulomb parameter, p the re¬ 
duced pp mass, C w 0.5772156649 the Euler number, and 
n the exponent in Eq. O- The renormalization phase 
4’r{p) to be used in the actual calculations with finite 
screening radii R is not unique, since only the infinite 


R limit matters, but the converged results have to show 
independence of the chosen form of (I)r{p)- According to 
our investigations this is indeed so. The results presented 
in this paper are based on the partial-wave dependent 
form iTl^ of the renormalization factor for which we 
find the convergence with R to be slightly faster than for 

( 1123 ). 

We refer to Refs. [nia for a rigorous justification of 
the renormalization procedure of Eqs. m and 03 and 
proceed here to study the numerical convergence of our 
predictions with increasing screening radius i? as a prac¬ 
tical justification for the validity of the chosen Coulomb 
treatment. 

The above discussion left out the identity of the two 
protons. Taking the identity of the protons into account, 
the transition amplitude {pfRf\t\piRi) of Eq. (Illbll has 
to be calculated for antisymmetrized states. Practical 
results based on Eq. (Illbll are shown in Figs. 121-El 

The Coulomb effect on the hadronic pp phase shifts p 
is most important in the ^S'o partial wave. The conver¬ 
gence with R for the ^S'o phase shift, shown in Fig. |21 
is impressive. The convergence is faster at higher ener¬ 
gies. A screening radius of i? = 20 fm (10 fm) suffices for 
an agreement within 0.01 deg with the exact phase shift 
values at all energies above 5 MeV (25 MeV). In con¬ 
trast, in order to reproduce the pp scattering length 
= —7.815 fm and the effective range = 2.773 fm 
within 0.010 fm, screening radii larger than R = 100 fm 
are required. In comparison to the screening function 
adopted in this paper, Fig.|21also proves the convergence 
with R to be rather slow for the Yukawa screening and to 
be of unpleasant oscillatory behavior for a sharp cutoff. 

Figure 121 studies the convergence of the result for the 
spin-averaged pp differential cross section at 5 MeV pro¬ 
ton lab energy with increasing screening radius R. The 
screening radius R = 20 fm appears to be sufficiently 
large for that energy, since, according to Fig. |21 the re¬ 
sults for i? > 20 fm are indistinguishable from the exact 
Coulomb results, despite the rather fine scale of the plot. 
The rate of convergence seen in Fig.|21is characteristic for 
all studied observables at that energy. The convergence 
of observables with R is also faster at higher energies; 
beyond 25 MeV the radius R = 10 fm is amply enough. 

Figure El shows the proton analyzing power results for 
pp scattering at 100 MeV proton lab energy. The re¬ 
sults are converged with respect to screening and the 
exact results are compared with two approximations, la¬ 
beled no Coulomb and Coulomb externally corrected: In 
the no Coulomb approximation, the Coulomb interaction 
is omitted completely; in the Coulomb externally cor¬ 
rected approximation, the Coulomb scattering amplitude 
is added to the no-Coulomb one, the latter being mod¬ 
ified for the external Coulomb distortion by multiplica¬ 
tion with the Coulomb phase factors in the ini¬ 

tial and final states M- Whereas the results converged 
with respect to screening and the exact results are indis¬ 
tinguishable in Fig. El the approximations no Coulomb 
and Coulomb externally corrected are pretty poor even 
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FIG. 2: Convergence of the pp phase shift p with screen¬ 
ing radius R for proton lab energies 3, 4, 5, and 10 MeV. Our 
results derived from Eq. lITT^ and given by dashed curves are 
compared with exact results given by solid lines. At 3 MeV 
also the results obtained with Yukawa screening (dashed- 
dotted curve) and with a sharp cutoff (dotted curve) are 
shown, demonstrating the superiority of the screening func¬ 
tion chosen in this paper. 
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FIG. 3: Convergence of the differential cross section for pp 
scattering at 5 MeV proton lab energy with screening radius 
R. The cross section is shown as function of the c.m. scatter¬ 
ing angle. Exact results given by the solid curve are compared 
to results with screening radius R = 10, 15 and 20 fm, given 
by dotted, dashed-dotted and dashed curves, respectively. Re¬ 
sults obtained with i? > 20 fm are not distinguishable from 
the exact results. 

at 100 MeV. Note that for the observable of Fig. O the 
results for both approximations lie out of the scale of that 
plot. 

The results presented in Figs. 13 and 0] are character¬ 
istic for all observables studied. We conclude that the 
employed method for the inclusion of the Coulomb inter¬ 
action in pp scattering works satisfactorily. We see con¬ 
vergence with increasing screening radius R at moderate 
values. The convergence in R is more rapid for higher 


FIG. 4: Proton analyzing power for pp scattering at 100 MeV 
proton lab energy as function of the c.m. scattering angle. 
Exact results, given by the solid curve, are indistinguishable 
from the results of our Coulomb treatment with screening ra¬ 
dius R > 10 fm; the dashed curve corresponds to no Coulomb 
results, while the dotted curve to the Coulomb externally cor¬ 
rected approximation. 

scattering energies; R = 10 fm is sufficient for proton lab 
energies above 25 MeV, whereas the screening radius is to 
be increased beyond 20 fm for energies below 5 MeV. We 
also note that the convergence in R is considerably slower 
for Yukawa screening and is of oscillatory behavior for a 
sharp cutoff. The exact Coulomb results are correctly ap¬ 
proached by the employed method with satisfactory ac¬ 
curacy, unlike the no Coulomb or the Coulomb externally 
corrected approximations. The method we use, based on 
the ideas of Refs. [Hill, encour^es us to carry it over 
to elastic pd scattering as Refs. P, Q did and to e.m. 
reactions involving the pd system. 

B. Elastic proton-deuteron scattering 

This section carries over the treatment of the Coulomb 
interaction, given in Sec. ESI for pp scattering, to elas¬ 
tic pd scattering. It establishes a theoretical procedure 
leading to a calculational scheme. 

Each pair of nucleons (/Jy) interacts through the strong 
coupled-channel potential Va and the Coulomb potential 
Wa- We assume that Wa acts formally between all pairs 
(/Sy) of particles, but it is nonzero only for states with 
two-charged baryons, i.e., pp and pA+ states. We intro¬ 
duce the full resolvent G^^\Z) for the auxiliary situation 
in which the Coulomb potential Wa is screened with a 
screening radius i?, Wa being replaced by WaR, 

G<^^\Z) = [Z-Ho-Y,v„-Y,w^r)-\ (13) 

<7 <7 

where Hq is the three-particle kinetic energy operator. 
The full resolvent yields the full pd scattering state 
when acting on the channel state |(()Q,(q)r'a) of relative 
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pd momentum q, energy Ea{q) and additional discrete 
quantum numbers Va and taking the appropriate limit 
Z = Ea{q) + iO. The full resolvent therefore also yields 
the desired S matrix. The full resolvent de¬ 

pends on the screening radius R for Coulomb and that 
dependence is notationally indicated; the same will be 
done for operators related to G^^^{Z). The full resol¬ 
vent G^^\Z), following standard AGS notation M of 
three-particle scattering, may be decomposed into chan¬ 
nel resolvents 

Gi^'>{Z) = {Z-Ho-Vc.-Wc.R)-\ (14) 

where, in pd channels a, WaR = 0, and into the full multi¬ 
channel three-particle transition matrix U^^{Z) accord¬ 
ing to 

G^^\Z) = 50^G\^HZ) + Gf\z)u'g{Z)G\^\Z). 

(15) 

The full multichannel transition matrix satisfies the AGS 
equation ^ 

= ho.G^\Z) + Y, h^Ti'^\z)Go{Z)Uil\z), 

(16a) 

where the two-particle transition matrix is derived from 
the full channel interaction Va + WaR, i.e., 

Ti^\Z) = (u„ + Wc,r) + (u„ + Wc,R)Go{Z)Tj^^'>{Z), 

(16b) 

Go{Z) = [Z — Ho)~^ being the free resolvent and dfSa = 
1 ~ S/ 3 a- Of course, the full multichannel transition ma¬ 
trix ujj^\Z) must contain the pure Coulomb transition 
matrix T^]^{Z) derived from the screened Coulomb po¬ 
tential between the spectator proton and the center 

of mass (c.m.) of the remaining neutron-proton {np) pair 
in channel a, 

{Z) = {Z)T^^^ (Z), (17) 

the pd channel being one of those channels a. This part 
may therefore be isolated in the full resolvent according 
to 

G(«)(Z) = Sp^Gi^Hz) + Gf\z)6p^T^^-iZ)Gi^Hz) 

+ Gjf^(Z)[t7^f)(Z) - 6p^T^^iZ)]Gi^Hz). 

(18) 

Nevertheless, as we have done in Sec. nm an alter¬ 
native decomposition of the full resolvent, which ap¬ 
pears conceptually neater for the purpose of elastic pd 
scattering, may be developed based on the following 
idea. Instead of correlating the plane-wave channel state 
l'('a(q)^'a) in a single step to the full scattering state by 
G*-^^ (Z), it may be correlated first to a screened Coulomb 
state of proton and deuteron by the screened Coulomb 


potential between a proton and the c.m. of an np 

pair through 

GMZ) = {Z-Ho-Vo.- w^R - (19a) 

G^r{Z) = Gi«)(Z) + Gi«)(Z)T=“(Z)Gi«)(Z), (19b) 

where, in each channel a, w^r and are never si¬ 

multaneously present: When a corresponds to a pp pair, 
WaR is present and = 0; when a denotes an np pair, 

WaR = 0 and is present. Thus, the full resolvent 

can alternatively be decomposed into 

G^^\Z) = 5p^G^r{Z) + G^r{Z)U^^I\z)G^r{Z), 

(20a) 

G(«)(Z) = S^^Gi^\z) + Gf\z)S^^T^SiZ)Gi^\z) 

+ Gf\z){[l + T^^t{Z)Gf\z)]uj,^\z) 

X [1 + Gi^\Z)T^liZ)]}Gi^Hz), (20b) 

where the operator UpJ{Z) may be calculated through 
the integral equation 

GpaiZ) = SpoXGakiZ) + Va] + Sfjoy^aR 

+ ^(^/3o•^'cr + 5i3aWl3R)GcrR{Z)U^^{Z), 

( 21 ) 

which is driven by the strong potential Va and the po¬ 
tential of three-nucleon nature WaR = J^ai^aaWaR — 
This potential WaR accounts for the differ¬ 
ence between the direct pp Coulomb interaction and the 
one that takes place between the proton and the c.m. of 
the remaining bound as well as unbound np pair. When 
calculated between on-shell screened pd Coulomb states, 
U^a\Z) is of short-range, even in the infinite R limit. 
Equation (I20bll gives an alternative form for the differ¬ 
ence of the transition matrices [uj^^ {Z) — Sj^aT^R (Z)] 
in Eq. (tUHl . i.e., 

Ui,^\z) - SpaT-^{Z) = [1 + T^^tiZ)Gf\z)]uj,^\z) 

X [l + Gi^Hz)T^^{Z)]. 

( 22 ) 

Though we calculate that difference directly from the 
potentials Uq, WaR, and through the numeri¬ 

cal solution of Eqs. m and ini), Eq. ® demon¬ 
strates that for initial and final pd states [G^^^(Z) — 
(Z)] is a short-range operator due to the nature 
of ujja'^iZ) as discussed above, but it is externally dis¬ 
torted due to the screened Coulomb wave generated by 
[1 -I- G^\Z)T^^(Z)]. Thus, Eq. ll^ achieves a clean 
separation of the full on-shell transition matrix (Z) 
into the long-range part dpaT^R'iZ) and the short-range 
part {Z) — 5i3aT^R{,Z)]. On-shell, both parts do not 
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have a proper limit as i? —> oo. They have to get renor¬ 
malized as the corresponding amplitudes for pp scattering 
in Sec. ESI in order to obtain the results appropriate for 
the unscreened Coulomb limit. 

According to Refs. P P, the full pd transition 
amplitude for initial and final states and 

|</'/3(q/)j^/3/); df = di, referring to the strong potential 
Va and the unscreened Coulomb potential Wa, is obtained 
via the renormalization of the on-shell multichannel tran- 
sition matrix (Z) with Z = Ea{di)+iO in the infinite 
R limit 

(0/3(q/)j^/3/|t^/3aka(q»);^ai) 

K —>^00 

X U^^^\Ea{d^) + iO)\<pa{q^)iyai)ZR^ id^)}■ 

(23a) 

As for pp scattering, the split of the full on-shell 
multichannel transition matrix into long- 

and short-range parts is most convenient. For the 
screened Coulomb transition matrix T^^{Z)^ contained 
in uj^^{Z), the limit in Eq. (I23all can be carried out 
analytically, yielding the proper Coulomb transition am¬ 
plitude {(l)f3{qf)iyi3f\T^c\M'ii)’^ai) PP, i-e., 

(0/3(Q/)^/3/ \ U0oi\4^a. 

+ " («/ ) ^^0 (q/ )^/3/ I {d^) + iO) 

- (Aa((7i) +*0)]|(/)„(q*)i^„,)Z”"(g,)}. 

(23b) 

The limit for the remaining part {Z) —5i3aT^^{Z)] 
of the multichannel transition matrix has to be performed 
numerically, but, due to the short-range nature of that 
part, it is reached with sufficient accuracy at finite screen¬ 
ing radii i?, and furthermore, {Z) — 5paT^R{Z)] can 

be calculated using a partial-wave expansion. 

In close analogy with pp scattering, the renormaliza¬ 
tion factor for i? ^ 00 is a diverging phase factor 

2fl(<z) (24a) 

where <^fi(g), though independent of the pd relative an¬ 
gular momentum I in the infinite R limit, is 

(l)R{d) ^ cri{q) - T]iR{q), (24b) 

with the diverging screened Coulomb pd phase shift 
ijiuid) corresponding to standard boundary conditions 
and the proper Coulomb one cri(q) referring to the log¬ 
arithmically distorted proper Coulomb boundary condi¬ 
tions. In analogy to pp scattering the form (I24bll of the 
renormalization phase is readily understood by looking 
back to Eq. (EHi. For the screened Coulomb potential of 
Eq. (E the infinite R limit of 4>Rid) is known analytically 

^Rid) = K(9)[ln (2gR) - C/n], (24c) 


K{q) = aM/q being the pd Coulomb parameter and M 
the reduced pd mass. The form of the renormalization 
phase Ip Rid) to be used in the actual calculations with 
finite screening radii R is not unique, but, like in Sec. EH 
the converged results show independence of the chosen 
form of (pRid)- The results presented in this paper are 
based on the partial-wave dependent form Ij24hll of the 
renormalization factor for which we find the convergence 
with R to be slightly faster than for 

Again we refer to Refs. PP for a rigorous justifi¬ 
cation of the correction procedure of Eqs. IIS and lEH) 
and proceed here to study the numerical convergence of 
our predictions with increasing R. 

We choose an isospin description for the three baryons 
involved in three-nucleon scattering. In the isospin for¬ 
malism the isospin T of the interacting pair and the 
isospin t of the spectator are coupled to the total isospin 
T with the projection A4t- Due to the hadronic charge 
dependence together with the screened Coulomb interac¬ 
tion in pp and pA'^ pair states, i.e., in states with isospin 
\TMt) = |11), the two-baryon transition matrix T^\Z) 
becomes an operator coupling total isospin T = ^ and 
T = I states according to 

{{T't')T' (Z) \ (TtjTAdr) 

= ST’TdftdM’^Mr ^ {TMTtmt\T'M t) 

MT'lTlt 

X {TMT\Ti^\Z)\TMT){TMTtmt\TMT)- 

(25) 

Due to the isospin formulation, the nucleons are therefore 
considered identical. However, the discussion has left out 
the identity of nucleons till now. Instead of the transi¬ 
tion amplitude of Eq. (I23bll we therefore have to use the 
properly symmetrized form 

iPa (q/ )^af I U\pcx iQ.i)Rai ) 

= '^{Pai(lf)Vaf\Uaa\Paiqi)l^cyi), 

<7 

(Q/)^a/ I^|0Q ) 

+ {ZPi^ idf){Pai^f)l'af\[U^^HEcid^) + iO) 

H —^oo 

- T^RiEaidi) + *0)]|«!)a(q*)j^c.,)2:”= (g,)} 

(26b) 

with C/(«)(Z) = uL^\z) + uifiZ)P23i + ui^\z)P 3,2 
for the calculation of observables, (a/dj) being cyclic and 
P 231 and P 312 being the two cyclic permutations of (a/dj). 
U^^'^iZ) satisfies the standard symmetrized form of the 
integral equation llball . i.e., 

U^^\Z) = PGQ^iZ) + PTi^\Z)Go{Z)U‘'^\z) (27) 
with P = P 231 + P 312 - 


The practical implementation of the outlined calcu- 
lational scheme faces a technical difficulty. We solve 
Eq. 1127|l in a partial-wave basis. The partial-wave expan¬ 
sion of the screened Coulomb potential converges rather 
slowly. The problem does not occur in pp scattering, 
since there the partial waves with different two-baryon 
total angular momentum I are not coupled and the max¬ 
imal I required for — tij(z)] is determined ac¬ 

cording to Eqs. and m by the range of the hadronic 
potential v. However, in the calculation of all 

two-baryon partial waves are coupled dynamically; the 
required maximal I is determined by the range of the 
screened Coulomb potential and is considerably higher 
than required for the hadronic potential alone. In this 
context, the perturbation theory for higher two-baryon 
partial waves developed in Ref. [13 is a very efficient and 
reliable technical tool for treating the screened Coulomb 
interaction in high partial waves. Furthermore, in prac¬ 
tical calculations we split the difference of the transition 
matrices in Eq. Il26bll into two parts with different partial- 
wave convergence properties, 

[/(«) {Z) - (Z) = (Z) - PT^r{Z)P] 

-[T^liZ)-PT^RiZ)P], 

TaR{Z) being the two-baryon screened Coulomb tran¬ 
sition matrix derived from WaR alone and hidden in 
Ti^\z) according to Eq. (llBbll . The term PTaR{Z)P 
is the remainder of the three-body operator U'^^\Z) — 
PGq^{Z) in the absence of the strong force, and it is 
contained in U^^\Z) as the most important Coulomb 
contribution; the difference \U^^\Z) — PTaR{Z)P] con¬ 
verges with respect to included two-baryon states consid¬ 
erably faster than U^^^{Z) alone. The term \T^'^-{Z) — 
PTaR{Z)P] accounts for the off c.m. pd screened 
Coulomb interaction and converges rather slowly, but 
the inclusion of very high partial waves is much easier 
than for [U^^'^{Z) — PTaR{Z)P]. We vary the dividing 
line between partial waves included exactly and pertur- 
batively in 1/(^1 [Z) as well as angular momentum cut¬ 
offs for both terms in Eq. (I28II in order to test the con¬ 
vergence and thereby establish the validity of the proce¬ 
dure. The problem of high partial waves does not occur in 
Refs. [10 , since the authors use the quasiparticle formal¬ 
ism and work with equations of two-body type in which 
the partial-wave decomposition has to be performed only 
with respect to the relative motion of the spectator par¬ 
ticle and the correlated pair. Due to technical limita¬ 
tions, Refs. dll use low-rank separable potentials for 
the hadronic interaction and approximate the two-proton 
screened Coulomb transition matrix by the potential. In 
contrast, we work with a realistic hadronic interaction 
(CD Bonn or CD Bonn -|- A) without separable approxi¬ 
mation, and we never approximate the energy-dependent 
pair transition matrix for screened Coulomb TaR{Z) by 
the potential WaR. 

With respect to the partial-wave expansion in the ac¬ 
tual calculations of this paper, we obtain fully converged 
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FIG. 5: Convergence of the differential cross section and of 
the proton analyzing power Ay{N) for pd elastic scattering 
at 3 MeV proton lab energy with screening radius R. The 
observables are shown as functions of the c.m. scattering 
angle. The hadronic potential is CD Bonn -|- A. Results ob¬ 
tained with screening radius i? = 5 fm (dotted curves), 10 fm 
(dashed-double-dotted curves), 15 fm (dashed-dotted curves), 
20 fm (double-dashed-dotted curves), 25 fm (solid curves) 
are compared. Results without Coulomb (dashed curves) are 
given as reference for the size of the Coulomb effect. 


results by taking into account the screened Coulomb in¬ 
teraction in two-baryon partial waves with pair orbital 
angular momentum L < 13 for the first term in Eq. |2H|) 
and with L < 25 for the second term; orbital angular mo¬ 
menta L > 7 can safely be treated perturbatively. The 
above values refer to the screening radius R = 25 fm; for 
smaller screening radii the convergence in orbital angular 
momentum is faster. The hadronic interaction is taken 
into account in two-baryon partial waves with total an¬ 
gular momentum / < 5. Both three-baryon total isospin 
T = i and T = | states are included. The maximal 

three-baryon total angular momentum J' considered is 
31 

2 ■ 

Figures [S] and 0 study the convergence of our 
method with increasing screening radius R according to 
Eq. The comparison with the no Coulomb results 

(dashed curves), used till now by us when accounting for 
pd data, gives the size of the Coulomb effect. First 
we concentrate on 3 MeV proton lab energy, the lowest 
energy considered in this paper. As examples we show 
the differential cross section, the nucleon analyzing power 
Ay{N) which has the most critical convergence behavior 
according to Refs. 0111 , and the deuteron tensor ana¬ 
lyzing power T 2 \, the most slowly converging observable 
at 3 MeV proton lab energy according to our experience. 
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FIG. 6: Convergence of the deuteron tensor analyzing power 
T 21 for pd elastic scattering at 3 MeV and 10 MeV proton 
lab energy with screening radius R. The observable is shown 
as function of the c.m. scattering angle. The curves are ex¬ 
plained in the caption of Fig. 

Nevertheless, the convergence is impressive even for those 
worst cases: Only T 21 , shown in Fig.|^ requires a screen¬ 
ing radius R > 15 fm. Convergence is more rapid at 
higher energies as demonstrated in Fig.|^for the deuteron 
tensor analyzing power T 21 at 3 and 10 MeV proton 
lab energy. The observed convergence strongly suggests 
the reliability of the chosen Coulomb treatment. Fur¬ 
thermore, the forthcoming paper M makes a detailed 
comparison between the results obtained by the present 
technique and those of Ref. Q obtained from the varia¬ 
tional solution of the three-nucleon Schrodinger equation 
in conhguration space with the inclusion of an unscreened 
Coulomb potential between the protons and imposing 
the proper Coulomb boundary conditions explicitly. The 
agreement, across the board, between the results derived 
from two entirely different methods, clearly indicates that 
both techniques for including the Coulomb interaction 
are reliable; this is another justihcation for the technique 
used in this paper. 

As in Fig.2]for pp scattering, Fig. [7| compares predic¬ 
tions including the Coulomb interaction with results from 
traditional approximate treatments which were labeled 
before as no Coulomb and Coulomb externally eorrected. 
As already known from Refs. UM both approximations 
are unsatisfactory at low energies. At higher energies the 
Coulomb effect is confined more and more to the forward 
direction; the no Coulomb treatment fails there, whereas 
the Coulomb externally corrected approximation is usu¬ 
ally not accurate enough for larger scattering angles. 

The seen Coulomb effects and their physics implica- 
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FIG. 7: Differential cross section and deuteron analyzing 
power Ay(d) for pd elastic scattering at 135 MeV proton lab 
energy as functions of the c.m. scattering angle. Converged 
results of the present Coulomb treatment with A = 10 fm 
given by solid curves are compared to the results calculated 
with no Coulomb (dashed curves) and Coulomb externally 
corrected (dotted curves) approximations. 

tions are discussed in Sec. nni 


C. Radiative capture and two-body e.m. 
disintegration of the three-nucleon bound state 

For the description of the considered e.m. processes the 
matrix element {ipa ^(q/)^'ay: K+)|R) of the e.m. 

current operator between the three-nucleon bound state 
and the pd scattering state has to be calculated. The 
calculation of that matrix element without Coulomb is 
discussed in great length in Refs. |^. This sub¬ 
section only discusses the modification which arises due 
to the inclusion of the Coulomb interaction between the 
charged baryons. Coulomb is included as a screened po¬ 
tential and the dependence of the bound and scattering 
states, i.e., |i?(^() and on the screen¬ 

ing radius R is notationally made explicit. In analogy to 
pd scattering, the current matrix element referring to the 
unscreened Coulomb potential is obtained via renormal¬ 
ization of the matrix element referring to the screened 
Coulomb potential in the infinite R limit 

(iAi-)(q/)i.„,lj^(Q,K+)|R) 

( 29 ) 
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FIG. 8: Convergence of the differential cross section and of 
the deuteron analyzing power r 2 i for pd radiative capture 
at 3 MeV proton lab energy with screening radius R. The 
observables are shown as functions of the c.m. scattering 
angle. The curves are explained in the caption of Fig.|^ 


As for pd scattering, the practical results presented in 
this paper are based on the partial-wave dependent form 
of the renormalization factor 124bll . Due to the short- 
range nature of j^(Q,K+)|i?(^() the limit i? ^ oo is 
reached with sufficient accuracy at finite screening radii 
R. The presence of the bound-state wave function in the 
matrix element strongly suppresses the contribution of 
the screened Coulomb interaction in high partial waves, 
i.e., two-baryon partial waves with orbital angular mo¬ 
mentum L < 6 are sufficient for convergence. The 
other quantum-number related cutoffs in the partial-wave 
dependence of the matrix element are the same as in 
Refs. PEI, i.e., I < 4, J' < ^ for photoreactions, 
and I < 3, J < for two-body electrodisintegration of 
^He. All calculations include both total isospin T = i 
and T = ^ states. 

Figure IHl studies the convergence of our method with 
increasing screening radius R for pd radiative capture 
at 3 MeV proton lab energy. We show the differential 
cross section and the deuteron tensor analyzing power 
T 21 which are the most critical observables in terms of 
convergence behavior. As in the case of pp and pd elas¬ 
tic scattering the convergence is impressive and becomes 
more rapid with increasing energy; it is quite comparable 
to pd elastic scattering. The convergence with increasing 
screening radius R is the same for two-body electrodis¬ 
integration of ^He; we therefore omit a corresponding 
figure. 


D. Conclusions on the practical implementation of 
the Coulomb interaction 

Using the described method we are able to include 
the Coulomb interaction between two protons in the de¬ 
scription of hadronic and e.m. three-nucleon reactions 
in the pd c.m. energy regime from about 1 MeV up to 
the pion production threshold. The screening radius re¬ 
quired for the convergence decreases with increasing en¬ 
ergy. Whereas i? = 20 fm is required for energies around 
deuteron breakup threshold, the screening radius can be 
lowered to i? = 10 fm above 10 MeV c.m. energy. In 
contrast, the screening radius has to be increased consid¬ 
erably when calculating extreme low-energy quantities, 
such as the pd doublet scattering length, which at present 
is outside the reach of our adopted technique. On the 
other hand, the high-energy limit is imposed by the form 
of the hadronic interaction which is applicable only be¬ 
low pion production threshold. We notice no particular 
feature of the convergence when crossing the three-body 
breakup threshold. However, the paper does not treat 
the three-body breakup reactions yet. 

III. RESULTS 

We base our calculations on the two-baryon coupled- 
channel potential CD Bonn -|- A with and without 
Coulomb and use the CD Bonn potential with and with¬ 
out Coulomb as purely nucleonic reference. We use the 
charge and current operators of Refs. [HEl , appro¬ 
priate for the underlying dynamics. We add relativistic 
corrections to the charge in the Siegert part of the oper¬ 
ator when describing the photoreactions, an admittedly 
questionable procedure, but entirely unrelated to the real 
issue of Coulomb in this paper. 

Obviously, we have much more predictions than it is 
possible and wise to show. Therefore we make a judicious 
selection and present those predictions which illustrate 
the message we believe the results tell us. The readers, 
dissatisfied with our choice, are welcome to obtain the 
results for their favorite data from us. 

Our predictions are dominantly based on the two- 
baryon coupled-channel potential CD Bonn -|- A; its 
single virtual A-isobar excitation yields, in the three- 
nucleon system, an effective three-nucleon force consis¬ 
tent with the two-nucleon interaction. A-isobar effects 
increase the ^H binding energy from 8.004 MeV for CD 
Bonn to 8.297 MeV for CD Bonn -|- A, the experimental 
value being 8.482 MeV. That binding energy increase has 
simultaneous beneficial effects on other bound state prop¬ 
erties, e.g., on the charge radius, but those effects also 
appear in the pd elastic scattering amplitude, especially 
in the three-nucleon = i"*" partial wave. The corre¬ 

lation between trinucleon binding and other low-energy 
observables is known as scaling. However, the Coulomb 
interaction also makes a significant contribution to trin¬ 
ucleon binding; the ^H — ^He binding energy difference is 
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0.746 MeV for CD Bonn and 0.756 MeV for CD Bonn + 
A, compared with the experimentally required value of 
0.764 MeV. This binding energy difference is dominantly 
due to the Coulomb repulsion between the protons in 
^He; the contribution arising from the hadronic charge 
asymmetry is much smaller. 

In three-nucleon scattering the A isobar therefore con¬ 
tributes to the scaling phenomenon at low energies, but 
it manifests itself more directly at higher energies when 
channel coupling becomes more probable. This sec¬ 
tion tries to explore the interplay between A-isobar and 
Coulomb effects in the considered three-nucleon reac¬ 
tions. 


A. Elastic proton-deuteron scattering 

Figures El and El give characteristic low-energy results. 
As examples we show observables at 3 MeV and 9 MeV 
proton lab energy, respectively below and above deuteron 
breakup threshold. The Coulomb effect is quite appre¬ 
ciable at all scattering angles, but its relative importance 
decreases with increasing energy. In contrast, on the 
scale of the observed Coulomb effect, the A-isobar ef¬ 
fect is minute at those low energies. The inclusion of 
Coulomb is essential for a successful account of data for 
the spin-averaged differential cross section and for the 
deuteron tensor analyzing powers. However, the inclu¬ 
sion of Coulomb increases the discrepancy between the¬ 
oretical predictions and experimental data in the peak 
region of proton and deuteron vector analyzing powers, 
the so-called A^-puzzle. Our findings are consistent with 
the results of Refs. HQ 

FigurelTTIshows selected results at 135 MeV proton lab 
energy. The Coulomb effect is confined to the forward 
direction, i.e., to c.m. scattering angles smaller than 30 
degrees where the A-isobar effect is not visible. The A- 
isobar effect shows up rather strongly in the region of the 
diffraction minimum, where its effect is beneficial and the 
Coulomb effect is gone. A-isobar and Coulomb effects are 
nicely separated. Thus, the A-isobar effect found previ¬ 
ously [13 on the Sagara discrepancy and on spin observ¬ 
ables remains essentially unchanged by the inclusion of 
the Coulomb interaction. The predictions of Fig. ITTlare 
characteristic for all observables at higher energies. 


B. Proton-deuteron radiative capture 

The e.m. current operator is the standard choice of 
Ref. [23 supplemented by the relativistic one-nucleon 
charge corrections, also given in Ref. |^, which we found 
to be important for some spin observables even at low en¬ 
ergies. 

References [23, [^ carried out corresponding realis¬ 
tic calculations for pd radiative capture with different 
two-nucleon potentials and an irreducible three-nucleon 


force, but without relativistic one-nucleon charge correc¬ 
tions. The calculations of Ref. take the Coulomb 
interaction fully into account, but are limited to reac¬ 
tions below 10 MeV c.m. energy. Reference [^ ne¬ 
glects the Coulomb interaction in the continuum states. 
When comparable, the results of this paper and those of 
Refs. [^[^ agree qualitatively; benchmark comparisons 
have not been done. 

Figure [El shows the Coulomb effect for pd radiative 
capture at 3 MeV proton lab energy. The Coulomb ef¬ 
fect is most important for the differential cross section 
which is reduced by about 20% and agrees rather well 
with the experimental data. In contrast, the spin ob¬ 
servables show only a small Coulomb effect. The effect of 
relativistic one-nucleon charge corrections is entirely neg¬ 
ligible for the differential cross section, but rather sizable 
and necessary for a satisfactory description of the data 
for the vector analyzing powers. Our results without rela¬ 
tivistic one-nucleon charge corrections are consistent with 
the corresponding calculations of Ref. [23 which also fail 
to account for the vector analyzing power data. A mod¬ 
erate A-isobar effect due to scaling is visible around the 
peak of the differential cross section. 

Selected deuteron analyzing powers at 17.5 MeV 
deuteron lab energy with moderate Coulomb effects are 
shown in Fig. together with the experimental data. 
Since the deuteron analyzing power Ayy is rather flat be¬ 
tween 40 deg and 140 deg according to Fig. El Fig- 
focuses on the energy dependence of Ayy at 90 deg pho¬ 
ton lab scattering angle. Clearly, our calculation ac¬ 
counts rather well for the known data of Ayy in the whole 
deuteron lab energy region up to 95 MeV. However, a 
similar study for T 20 at 90 deg (not shown here) indicates 
that the strong energy dependence of the low energy data 
from TUNL is not compatible with the present cal¬ 
culation. The rather good agreement with experimental 
data in Figs. El and El is obtained in general as an in¬ 
terplay of three considered effects, i.e., the effects due to 
the A isobar, due to the relativistic one-nucleon charge 
corrections, and due to the Coulomb interaction. 

Figure El shows the differential cross section and the 
nucleon analyzing power for pd radiative capture at 
150 MeV nucleon lab energy where we previously found 
rather significant A-isobar effects [^. Even at this 
relatively high energy there is a visible, though small 
Coulomb effect around the peak of the differential cross 
section. In addition, both observables show a sizable ef¬ 
fect of the relativistic one-nucleon charge corrections. 

Figures El- cni also recall the non-Coulomb effects 
on observables in order to discuss their interplay with 
Coulomb. In the differential cross section at 3 MeV and 
150 MeV proton lab energy the A isobar plays different 
roles; at low energy the A isobar manifests itself through 
scaling due to changed bound state properties, but at 
higher energies the explicit excitation of A channels in 
scattering becomes more predominant. We also note 
that relativistic one-nucleon charge corrections are im¬ 
portant for proton and deuteron vector analyzing powers 
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FIG. 9: Differential cross section and analyzing powers for pd elastic scattering at 3 MeV proton lab energy as functions of 
the c.m. scattering angle. Results including A-isobar excitation and the Coulomb interaction (solid curves) are compared to 
results without Coulomb (dashed curves). In order to appreciate the size of the A-isobar effect the purely nucleonic results 
including Coulomb are also shown (dotted curves). The experimental data are from Ref. 1^ . 








FIG. 10: Differential cross section and analyzing powers for pd elastic scattering at 9 MeV proton lab energy as functions of 
the c.m. scattering angle. The curves are explained in the caption of Fig. M The experimental data are from Ref. 1^ . 
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FIG. 11: Differential cross section and proton analyzing 
power Ay (N) for pd elastic scattering at 135 MeV proton lab 
energy as functions of the c.m. scattering angle. The curves 
are explained in the caption of Fig.|^ The experimental data 
are from Ref. (crosses) and from Ref. (full circles) 
for the differential cross section, and from Ref. 1^ for the 
analyzing power. 


even at low energies; with increasing energy they become 
quite significant in general. The relativistic one-nucleon 
charge corrections are mostly beneficial for the account of 
the experimental data, though their inclusion is not fully 
consistent with the underlying nonrelativistic hadronic 
dynamics. 


C. Two-body electrodisintegration of ®He 

The e.m. current operator is taken over from Ref. |2l|: 
compared to photo reactions, the relativistic one-nucleon 
charge corrections are less important and are therefore 
omitted. The Coulomb effect in the two-body electro¬ 
disintegration of ^He depends on both energy and mo¬ 
mentum transfer. We do not study that dependence in 
detail. We only show in Fig. a sample result for the 
three reaction kinematics Cl, I and HR of Ref. The 
Coulomb effect on the Cl and I differential cross sections 
is visible in the peak, though in Cl small compared with 
the discrepancy between theoretical predictions and ex¬ 
perimental data. In the HR differential cross section a 
Coulomb effect is not visible in the logarithmic scale of 
the plot, but instead a A-isobar effect at backward an¬ 
gles. Qualitatively our results without Coulomb agree 
well with the ones of Ref. . 


IV. SUMMARY 

This paper shows how the Coulomb interaction be¬ 
tween the charged baryons can be included into the 
momentum-space description of elastic proton-deuteron 
scattering and of related e.m. reactions using the 
screening and renormalization approach. The theoreti¬ 
cal framework is the AGS integral equation M The 
calculations are done on the same level of accuracy and 
sophistication as for the corresponding neutron-deuteron 
reactions. The conclusions of the paper refer to the devel¬ 
oped technique and to the physics results obtained with 
that technique. 

Technically, the idea of screening and renormalization 
is the one of Refs. Pin and we rely on these works for 
mathematical rigor. However, our practical realization 
differs quite significantly from the one of Refs. HQ: 

(1) We use modern hadronic interactions, CD Bonn 
and CD Bonn -f AHn contrast to the low-rank separable 
potentials of Refs. HQ. Our use of the full potential re¬ 
quires the standard form of the three-particle equations, 
different from the quasiparticle approach of Refs. HQ 

(2) We do not approximate the screened Coulomb tran¬ 
sition matrix by the screened Coulomb potential. 

(3) The quasiparticle approach of Refs. HQ treats the 
screened Coulomb potential between the protons with¬ 
out partial-wave expansion and therefore has no prob¬ 
lems with the slow convergence of that expansion. Our 
solution of three-nucleon equations proceeds in partial- 
wave basis and therefore is faced with the slow partial- 
wave convergence of the Coulomb interaction between 
the charged baryons. However, we are able to obtain 
fully converged results by choosing a special form of the 
screening function and by using the perturbation theory 
of Ref. for treating the screened Coulomb transition 
matrix in high partial waves. This would not be possible 
if we had used Yukawa screening as in Refs. HQ- 

(4) Our method for including the Coulomb interaction 
is efficient. Though the number of the isospin triplet 
partial waves to be taken into account is considerably 
higher than in the case without Coulomb, the required 
computing time increases only by a factor of 2 to 3, due to 
the use of the perturbation theory for high partial waves. 

The obtained results are stable and well checked for 
their validity. The employed technique gets cumbersome 
when approaching very low energies, i.e., pd c.m. energies 
below 1 MeV, due to the need for very large screening 
radii. Thus, the calculation of the doublet scattering 
length for elastic pd scattering is, at present, outside of 
our numerical reach, a barrier which does not exist for the 
coordinate-space techniques adopted in Refs. QQ On 
the other hand, we do not see any particular numerical 
problem when crossing the breakup threshold and going 
to higher energies where coordinate-space techniques are 
very hard to apply. 

The present technique is not yet used for breakup itself, 
but such an extension is on its way. 

Physicswise, the Coulomb effect in elastic pd scatter- 
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FIG. 12: Differential cross section and analyzing powers for pd radiative capture at 3 MeV proton lab energy as functions 
of the photon c.m. scattering angle with respect to the direction of the proton. The shown results include, respectively: 
Coulomb interaction, relativistic one-nucleon charge correction, and A-isobar excitation (solid curves); relativistic one-nucleon 
charge corrections and A-isobar excitation (dashed-dotted curves); A-isobar excitation (dashed curves); purely nucleonic results 
(dotted curves). The experimental data are from Ref. la^. 






FIG. 13: Deuteron analyzing powers for pd radiative capture at 17.5 MeV deuteron lab energy as functions of the photon 
c.m. scattering angle with respect to the direction of the proton. The curves are explained in the caption of Fig. \n\ The 
experimental data are from Ref. P. 
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FIG. 14: Deuteron analyzing power Ayy for pd radiative 
capture at 90 deg photon lab scattering angle as functions 
of the deuteron lab energy. The curves are explained in the 
cap tion of Fig. EH The experimental data are from Refs. 

il. 



0 60 120 180 


©c.m. (deg) 

FIG. 15: Differential cross section and proton analyzing 
power Ay{N) for pd radiative capture at 150 MeV proton 
lab energy as functions of the photon c.m. scattering angle 
with respect to the direction of the proton. The curves are 
explained in the caption of Fig. EH The experimental data 
are from Ref. lo^ . 


ing is important at low energies for all kinematic regimes, 
but gets confined to the forward direction at higher ener¬ 
gies, whereas the effect mediated by the A isobar remains 
almost unmodified by the inclusion of Coulomb. In ra¬ 
diative pd capture the Coulomb effect is important for 
low-energy differential cross sections and for some spin 




FIG. 16: Lab differential cross section for two-body electro¬ 
disintegration of ®He as function of the proton lab scattering 
angle. The electron lab energy, scattering angle, the momen¬ 
tum and energy transfer are 390 MeV, 74.4 deg, 434.8 MeV 
and 66.1 MeV for the reaction kinematics Cl, 527.9 MeV, 
52.2 deg, 430.0 MeV and 99.8 MeV for the reaction kinemat¬ 
ics I, and 390 MeV, 39.7 deg, 250.2 MeV and 113.0 MeV for 
the reaction kinematics HR of Ref. 1^ . respectively. Results 
including A-isobar excitation and the Coulomb interaction 
(solid curves) are compared to the results without Coulomb 
(dashed curves). In order to appreciate the size of the A- 
isobar effect the purely nucleonic results with Coulomb are 
also shown (dotted curves). The experimental data are from 
Ref. 1^. 
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observables up to about 30 MeV proton lab energy; at 
higher energies there is still a visible Coulomb effect for 
some observables, e.g., in the peak of the differential 
cross section. In two-body electrodisintegration of ^He 
the Coulomb effect appears not to be simply related to 
the internal excitation of the three-nucleon system. A 
thorough study of the dependence of the Coulomb effect 
on the energy- and three-momentum transfer to the ^He 
target is beyond the scope of this paper. 
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